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Abstract. A transport methodology to study the electron transport between 
quantum dots arrays based in Transfer Hamiltonian approach is presented. The 
interactions between the quantum dots and between the quantum dots and the 
electrodes are introduced by transition rates and capacitive couplings. The 
effects of the local potential are computed within the self-consistent field regime. 
The model has been developed and expressed in a matrix form in order to 
make it extendable to larger systems. Transport through several quantum dot 
configurations have been studied in order to validate the model. Despite the 
simplicity of the model, well-known effects are satisfactorily reproduced and 
explained. The results qualitatively agree with other results obtained using more 
complex theoretical approaches. 
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1. Introduction 

Confined structures have been available to the experimentalist for a very long time, 
the MOS (metal-oxide-semiconductor) transistor is the archetype of a confined two- 
dimensional system pQ. Nevertheless, the possibility to enhance this confinement 
by embedding low-dimensional structures in an insulating matrix has renewed the 
interest. These structures (quantum dots, wires or layers) can be used in single- 
electron device [2] ,new memory concepts |3j and photon or electroluminescent devices 

Concerning quantum dots (Qds), they are particularly attractive because they possess 
discrete energy levels and quantum properties similar to natural atoms or molecules. 
From a fundamental point of view, research has been mostly concentrated on single 
quantum dots. These simple systems have been studied using many-body approaches, 
including non-equilibrium Green's function formalism (NEGFF) |6]. From a 
practical point of view, many novel phenomena have been discovered, such as the 
staircaselike current- voltage (FV) characteristic [7], Coulomb blockade oscillation [5], 
negative differential capacitance [9J, and the Kondo effect in Qds [TU] , 
Researchers have recently paid much attention to electron transport through several 
Qds since multiple Qd provides more Feynman paths for the electron transmission 
However, up to now the only computation of transport in an extended arbitrary 
array of Qds was done by Carreras et al. [H] but no local potential due to self-charge 
was included. Sun et al.[T3] have also studied the electron transport using NEGFF for 
different arrangements of Qds, from one to three Qds, without including the potential 
due to self-charge neither. The inclusion of the self-charge potential using this complex 
framework is usually impossible for large systems. 

In this work, we use non-coherent rate equations (NCRE) [IU[T5] to study the electrical 
transport in Qds in an extendible, arbitrary, matrix of Qds taking into account self 
charge effects. In a previous work [TB], we applied NCRE to obtain analytical solutions 
for electron transport in simple cases. Using this approach each Qd is treated as a 
separate system, therefore we can write a NCRE for each dot since these equations 
describe relationship between the charge inside the Qds and applied bias voltage. The 
interactions between the Qds, and between the Qds and the electrodes are introduced 
by transition rates and capacitive couplings. Electron transport and charge densities 
inside the Qds depend on the tunnel transparency of the barriers limiting each dot. In 
order to effectively solve the multielectron problem, the effects of the local potential are 
computed within the self-consistent field (SCF) regime. Moreover, we show how our 
approach can be easily extended to an arbitrary number of Qds and configurations 
using a matrix formalism. Therefore, this methodology allows to simulate realistic 
devices based on large scale Qds arrays. Finally, we compare this methodology with 
NEGFF, obtaining similar results. 

2. Theoretical background 

Our system consist of two electrodes (L lead and R lead) coupled to a central transport 
region. The central region contains several quantum dots, N Qds, distributed inside of 
an insulator matrix. In order to find the current voltage curve I-V of the total system 
we use the transfer Hamiltonian formalism |171 ITS] . Using this formalism we can write 
an expression for the current flowing across two parts of the system. Assuming no 
inelastic scattering and symmetry in the transmission coefficient [TH] the net current 
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flux between two parts of the system is 
_ iwq 



TnWpiWpjiEKfjiE) - fi(E))dE, 



(1) 



where Tij (E) is the transmission probability, pi (E) and pj (E) are the density of states 
while fi (E) and fj (E) are the distribution functions of the different parts of the system. 
In the equilibrium, the electrochemical potential of the whole system is equal and the 
particular distribution functions (DFs) are described by the equilibrium Fermi Dirac 
DF therefore the current between each part of the system is zero. If an external 
bias voltage (V) is applied, which will drive the system out of the equilibrium, the 
electrochemical potential of the leads will change by pl—P-r. = QV. From the definition 
of the total charge Ni inside the i th Qd, we can write 



N 



Pi (E)ni(E)dE, 



(2) 



( / T Ll p L pi(f L - ni)dE + / T Ri p R pi(f R - ni)dE 



where rii(E) is an unknown DF and /?, is the density of states (DOS) of the i th Qd. 
For sake of clarity, we only consider one single state with energy level e in each Qd. 
In order to take into account the coupling with the surrounding elements we assign a 
Lorentzian shape DOS centered in e. We can write the evolution charge in time for 
each Qd as Ni — J Ijidt, where the subscript i refers to i th Qd and j runs over 
the other components of the system. Thus, a set of integro-differential equations are 
obtained for the time charge evolution 

dNi _ 4nq i 
~dT 
v 

+ ^2 TjiPjPiinj - Ui)dE) \ti = l...N, (3) 

where we explicitly write all the current terms: the leads current contributions (first 
and second term) and the neighbor contribution (the last term). We assume that the 
DFs in the electrodes (/r, and f R ) are similar to the Fermi Dirac DF using different 
electrochemical potentials (pl and pr). Equation ^ can be rewritten for the steady 
state and assuming no inelastic scattering we can obtain the DF in each Qd for each 
energy step as a solution of the system of equations 



/ -T L1 p L - T R1 p R - 



(N-l) 



T ijPj 



Tin Pn 



\ ( 



Tin Pi 
-ThiPhfL — T R ip R f R 



-T LN p L -T RN PR-J2j^N ) TnjPj ) \ n N 



(4) 



-ThNPhfh — T R n p R f R 

The effect of the applied voltage to the external electrodes on the electrostatic 
potential inside each Qd must also be taken into account. The classical solution for 
the potential at each quantum dot (Vi) involves the Poisson equation 

qANt 



rk 



(5) 



A Transfer Hamiltonian approach in self- consistent field regime for transport in arbitrary quantum dot arrays 4 



where e r is the relative permittivity of the dielectric media, eo is the vacuum 
permittivity and is the Qd volume. The general solution for the potential energy 
Ui = -qVi in the i th Qd is gD] 

where the subscript j runs over all components of the system, Cy is the capacitive cou- 
pling between the different components and C to t,i = J2j j^i ^ s the total capacitive 
coupling of i th Qd. The charge energy constant Uot — q 2 /Ctot,i is the potential increase 
as a consequence of the injection of one electron into the Qd and AJV; is the change 
in the number of electrons, calculated respect to the number of electrons N initially 
in the i th Qd. The effects of local potential on each Qd, which modify the Qd charge 
and the currents, should be taken into account in the Qd DOS Pi(E) —> p%(E — Ui). 
In (|6| we observe that the local potential depends on the increasing charge density 
but at the same time the charge depends on the DOS that it is modified by the local 
potential. These considerations impose a self-consistent solution of and (Rjl. 



3. Results and discussion 

In this section, we first show the calculated current voltage curves I-V for different 
arrangements and we compare their with the results obtained using NEGFF 13 . We 
also present the number of electrons JV^ accumulated in the i th Qd in each configuration. 
In these cases, analytical expressions for the current are presented as well. Finally, the 
extension of the model developed in the previous section is presented as a powerful 
method to study the electron transport in an arbitrary extended array of Qds. 
The electrochemical potentials in the two leads are set at pi = and pn = —qV. 
Electrons flow from the left lead to the right one. For simplicity, we consider that the 
transmission probability is constant and the same between all the parts of the system. 
We do not consider direct transmission between the leads. For clarity the DOS of 
the leads were considered constant in all the energy range. Using this framework, 
transport without inelastic scattering, the position of the energy levels in the Qds plays 
an important role therefore the evolution of it with the applied bias voltage define the 
shape of the I(V) curve. As expected, the I(V) curves exhibit strong dependence with 
the electrostatic coupling of the different parts of the system. We present expressions 
for the evolution of the energy level with the applied bias voltage assuming that the 
elements which are coupled have equal capacitive coupling between them. We set a 
constant charge energy for all Qds, Uq = 0.25eV. 

3.1. One single Qd 

We briefly review electron transport through one Qd. Using ([3| and taking into 
account only lead contributions the current can be written as 

I=~f- Tf, —7?, (fL-fR)dE. (7) 

n J J-LlPL + -LRlPr 

Fig. [I] shows the numerical result of the current I(V). In the calculation we assumed 
symmetric coupling respect to the leads, Tri = Tlx = 0.2 [|] The evolution of the 

X We use similar transmission values than Sun et al.|13] in order to make possible the qualitative 
comparison between models. 
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Figure 1. (a) The I-V curve for one single Qd obtained using NCRE. We also 
show the NEGFF results for the same system, the NEGFF data are taken from 
Sun et al. 113) - (b) The electron number in the Qd as a function of the applied 
bias V. The inset shows the connection geometry. The rectangles represent two 
leads and the circle represents a Qd. 



energy level with the applied bias voltage is 

e 1 (V) = l-V/2 + U AN 1 , (8) 

where the second and third terms are due to the electrostatic effect. As expected the 
current increases with the bias when the energy of the Qd moves across the left lead; 
which is fii = e\(V) — > V ~ 2. When V is bight enough, the current saturates to a 
constant value as £i(V) is placed between the two electrochemical potentials of the 
leads. Fig. U\b) shows the dependence of the electron number with the applied bias. 

3.2. Two Qds 

We now study the case of two Qds. There are four different connection geometries 
between Qds and leads. In our calculations we assume symmetric coupling respect to 
the leads, T m = T L1 = 0.2, and the Qd coupling T 12 = 0.2. 

3.2.1. Parallel case The first configuration of two Qds is the case that they are in 
parallel. Both Qd are coupled to all elements of the system, the leads and the neighbor 
Qd. In this configuration the expressions for the current are 

_ Airq f T L1 T R1 (T L1 p L + Ti R p R + T 12 (pi + p 2 ))phPRP\ 
h -^J D 2 

x(fL-f R )dE (9a) 
_ Airq f T L1 T R1 (T L1 p L + T 1R p R + T 12 (pi + P2))plPrP2 

l2 --fTJ D 2 

x(f L -f R )dE, (9b) 

where D 2 = (Ti R p L +T 1L p R ) 2 +T 1L T 12 p R (p 1 + p 2 ) +T L1 T 12 p L (p 1 +p 2 ). The position 
of the energy level of each Qd is 

e 1 (V) = 1 -qV/3-qV 2 /3 + U AN 1 (10a) 
e 2 (V) = 3.5 - qV/3 - qVi/3 + U Q AN 2 . (106) 
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Figure 2. (a) The total I-V curve and partial I-V curves obtained using NCRE 
for a parallel configuration. The NEGFF results are taken from Sun et al. 1 131 ■ 
(b) The electron number in the Qds as a function of the applied bias V. 



We show the total and partial currents Fig. [2]ja) . The I-V curve shows two steps when 
the energy levels of the Qds are placed between the electrochemical potentials of the 
leads. This case is equivalent to a single Qd with two energy levels. Fig. [2jb) shows 
the electron number ni with the applied bias voltage. The charge increases until it 
reach the saturation value. 
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Figure 3. (a) The I-V curve for two Qds in a serial configuration obtained using 
NCRE. We also show the NEGFF results for the same system, the NEGFF data 
are taken from Sun et al. 1 131 . The inset shows the connection geometry, (b) The 
electron number in the Qds as a function of the applied bias V. 



3.2.2. Serial case The second type of arrangement is the case of two Qds in a serial 
configuration. The system is shown in the inset of Fig. [3^a). Each Qd only interacts 
with one lead and the other Qd. In this case, the expression for the current is 

_ f T l1 Ti 2 T 2 rPlP\P2Pr 

h J T L iTi 2 pipL + T l1 T 2 b.PrPl + T 12 T 2R p 2 p R 
x(fL-f R )dE (11) 
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and the evolution of the energy level of each Qd with the applied bias voltage is 

e 1 (V) = l-qV 2 /2 + U AN 1 (12a) 

e 2 (V) = 3.5 - qV/2 - qVx/2 + U AN 2 , (126) 

where we assumed that the Qds are only coupled to each other and to one lead. In 
order to have current flowing through the system, the energy levels must lie between 
the electrochemical potentials of the leads and overlapping of the Qd energy levels 
is necessary. This means that the electrons need available states in each part of the 
system in order to move from Left lead to Right lead. When the energy levels are equal, 
£i = £2 V rs 7.5, this is a maximum overlapping between Qd DOS, the current is 
maximum and the system is in a resonance state therefore the channel is open. When 
the voltage increases further the Qd DOS overlapping decreases. Therefore a negative 
differential resistance appears [21] . In Fig. [3jb) we show the evolution of the charge 
in each Qd Nj as a function of the applied voltage V. Initially, N\ increases since the 
channel between the first and second Qd is closed. At the resonant condition, the 
channel between the Qds opens and some charge stored in the first Qd flows to the 
second Qd. At higher voltages the channel closes again and Ni stores all the incoming 
charge, while N 2 loses its charge. 
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Figure 4. (a) The I-V curve, for the configuration plotted in the inset, obtained 
using NCRE. We also show the NEGFF results for the same system, the NEGFF 
data are taken from Sun et al. 1131 . (b) The electron number in the Qds as a 
function of the applied bias V. 




3.2.3. Other two Qds configurations We first examine the case in which one Qd 
interacts with the two leads and it is also connected to the second Qd, while the 
second Qd is only connected to the first Qd. The current is 

n J -LliPl + -LRlPr 
and the position of the energy levels are 

ei(V) = 1 -qV/3-qV 2 /3 + U AN 1 (14a) 
e 2 (V) = 3.5 ~ qVi + U Q AN 2 . (146) 



The obtained current expression 13 is the same than the one we obtained for the single 
Qd case. The DF in the second Qd is the same as in the first Qd therefore the current 
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Figure 5. (a) The total I-V curve and partial I-V curves obtained using NCRE 
for the configuration showed in the inset. The NEGFF results are taken from Sun 
et al. I13| . (b) The electron number in the Qds as a function of the applied bias 
V. 



between the Qds is zero. The results are presented in Fig. |4j 

The second arrangement of Qds is shown in the inset of Fig. [5] The expressions 
for the current are 

Airq f Ti R pip R (T m p R T L1 p L +T 12 piT L1 p L ) t \jn(-\K\ 

h = J -f) v L ~ JR) dE y lba ! 

Airq [ T 2R p 2 p R T 12 piT L1 p L 
h = -j£- 7j (h - JR)dE, (156) 

where D = T 2R p R T L1 p L + T R1 T R2 p R + T R2 p R T 12 p 2 + T 12 p\T L1 p L + T R1 p R T 12 px and 
the total current is I = I\ + 1 2 . The energy level position is 

ei(V) = l-qV/3-qV 2 /3 + U AN 1 (16a) 
e 2 (V) = 3.5-qV 1 /2 + U AN 2 . (166) 

In this case we show the total and partial currents. The I-V partial current shows 
an interesting behavior. The current through the first Qd is similar to the single one 
Qd configuration but the current through the second reminds the slope of a resonant 
state. This fact can be easily understood in the following way: if the channel between 
the two Qds is closed the current only flows through the first Qd. When the Qdl-Qd2 
channel is opened the Qd2 also conducts. In the same case as before, when the voltage 
increases the overlapping decreases and the Qd2 current decreases. 



3.3. Three Qds 

The methodology developed in the first section can be easily extended into more 
complicated systems. Here, we present the results for some configurations based in 
three Qds. The analytical expressions for the current are too large to write here but in 
Fig.|6][(a),(c),(e)] we show the I-V curves and the charge in each Qd Fig.|6[(b),(d),(f)]. 
As we have shown before the position of the energy levels plays an important role in 
the I-V and N-V curves, using [6] we can write the position of the each energy level as 
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Figure 6. [(a),(c) and (e)] The I-V curves and the electron number [(b), (d) and 
(f)] respectively for three Qds with different configuration. The insets show the 
connection geometry. The NEGFF results are taken from Sun et al. |13| . 



a function of the applied bias voltage 

e 1 (V) = l-Y / -^-V j + U AN 1 (17a) 

3 

e 2 {V) = 2-Y j -^-V j + U /\N 2 (176) 

3 

es(V) = 3.5 - V p^-Vj + U AN 3 , (17c) 

^tota.13 

3 

where the subscript j runs over all connected elements of the system. The Qd-lead 
coupling and the interdot coupling are set equal Ty = 0.2. In the insets of the Fig. [6] 
we show the scheme of the system under study. 



3. 4- Large Qds arrangement 

To conclude we present the results for larger systems that they are close to the 
experimental measurements. The systems are formed by 100 Qds placed in a parallel 
configuration, serial configuration and in an array geometry (10 x 10). The total 
I-V curves and the geometries are presented in Fig. [7| The Qd-lead coupling and 
the interdot coupling are set equal Ty = 0.2. The capacitance between the linked 
elements are also equal. In order to represent an experimental system we considerer 
that the value of the energy level of each dot follows a normal distribution with mean 
value leV and deviation 0.2eV. This fact represents the usual distribution size that 
appears in the experiments. The relationship between the Qd radius and the energy 
level position is a well known effect and it is related to the quantum confinement of 
the electrons [22] . 

The I-V curves show an interesting behavior. First, in the parallel case, Fig. [7|a), the 
I-V curve shows a staircaselike structure and saturates to a constant value at high bias. 
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Figure 7. The I-V curves for the larger systems: (a) 100 Qds in parallel 
configuration, (b) 100 Qds in serial configuration and (c) 100 Qds in an array 
disposition 10 X 10. 

As we have seen before in the parallel configuration each Qd acts as an independent 
channel therefore the total current is the sum of all partial currents. As expected, the 
saturation current is 100 times the saturation current of a single dot. 
For the serial configuration, Fig.[7]Jb), we obtain a current peak as we expect due to the 
resonant state is necessary in order to have electron transport in this configuration. 
The maximum value of the peak is hard to determine because it depends on the 
transmission coefficient, but it also depends of the overlapping between the DOS of 
the Qds. 

Concerning to the array configuration, Fig. [7jc), the I-V curve is determined by a 
combination of the two previous cases. In order to have transport the resonant state 
condition must be fulfilled therefore a current peak appears but the total current is 
the sum of the partial currents of each row. 

3.5. Comparison with NEGFF 

Finally, the results obtained using the proposed approach have been compared to the 
results of Sun et al. |13j . In that paper the authors have used the nonequilibrium 
Green's function method (NEGFF) to study the electron transport between one, two 
and three Qds in several configurations. Their I-V results have been plotted in our 
figures (NEGFF in the legend). The main results are: 

• The results presented in Figs. |TJ [3j [2j |6ja) and |6jc) are in accordance between 
the two approaches. For the serial configuration Fig. [3] the differences are due 
to the different values of the Qd coupling, we also obtain a resonant peak when 
the energy levels of the Qds are placed in a resonant state. The resonant state is 
strongly dependent on the capacitive coupling of the Qds, as the position of the 
energy level with the applied bias voltage depends on the capacitive coupling of 
the Qd. 

In the parallel configuration we obtain the same staircase shape, but, in our case 
we also take into account the energy charge terms, therefore the steps occur at 
higher voltages. 

• The main difference appears in the case described in Fig. [4] For this configuration 
Sun et al. predicts an antiresonance effect. We do not recover this effect because 



A Transfer Hamiltonian approach in self- consistent field regime for transport in arbitrary quantum dot arrays 11 

our model considers each Qd as a separate quantum systems. For this reason, 
our approach is known as a non-coherent model. 

• For the systems presented in Figs. [5] and ^b) we obtain similar results. The 
position ol the current peak is different because Sun et al. assume that the bias is 
uniformly applied throughout the whole system meanwhile we take into account 
all the electrostatic coupling between the different parts of the system. 

As we have shown the electrostatic coupling plays an important role to determine 
the I-V characteristic of the system. The electrostatic effect has two terms: the first 
term is determined by the influence of the leads and the neighbor Qd and it is described 
by the capacitive coupling of the Qd and its surrounding. The second term takes into 
account the charge stored inside the Qd, this effect is related to the electron- electron 
interaction and the self consistent solution of[2]and[6]is the first approach to introduce 
many body effects, like the Coulomb Blockade. If we create nanodevices in oder to 
take advantage to the quantization of the current, only small number of discrete energy 
levels are available for conduction, the accurate control of the energy levels with the 
applied bias voltage is one of the most important points that we need to take into 
account. Therefore a good modelization of the Qd-Qd and Qd-lead capacitances is 
necessary. 

This paper precedes future works in which realistic DOS, energy dependent 
transmission coefficients as well as a realistic capacitive couplings can be introduced. 

4. Conclusion 

We propose a theoretical model to study the electron current in systems based in 
quantum dots (Qds). This model is based in the transfer Hamiltonian formalism and 
computes the I-V and N-V curves in the self consistent held regime (SCF), using non- 
coherent rate equations (NCRE). This approach provides a simple and transparent 
method to describe the electron transport. Due to the simplicity of the model, this can 
be easily extended to analyze arbitrary large arrays of Qds of interest in technological 
applications. Despite its simplicity and in contrast with other approaches the effect 
of self-charge has been taken into account, by solving the Poisson equation with 
appropriate boundary conditions for each Qd. As expected, the calculation of the 
local potential inside each Qd is one of the most critical points, since the I-V curves 
depends on the position of the energy level. 

In order to show the potential of this method to analyze realistic configurations, we 
have studied the electron transport between different Qd configurations. We have also 
compared the NCRE results with well established data obtained with the NEGFF 
approach. Such a successful comparison shows that NCRE is a powerful and intuitive 
method to describe the electron transport. 
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